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We study Josephson effect of d- wave superconductor (DS)/ferromagnet insulator(Fl)/DS junc¬ 
tions on a surface of topological insulator (Tl). We calculate Josephson current 1(f) for various 
orientations of the junctions where f is the macroscopic phase difference between two DSs. In certain 
configurations, we find anomalous current-phase relation 1(f) = — I (—f + n) with 27r periodicity. 

In the case where the first order Josephson coupling is absent without magnetization in FI, 1(f) 
can be proportional to cosy). The magnitude of the obtained Josephson current is enhanced due 
to the zero energy states on the edge of DS on TI. Even if we introduce an s-wave component of 
pair potential in DS, we can still expect the anomalous current-phase relation in asymmetric DS 
junctions with I (f = 0) ^ 0. This can be used to probe the induced d- wave component of pair 
potential on TI surface in high-T c cuperate/TI hybrid structures. 

PACS numbers: 74.45.+C, 74.50.+r, 74.20.Rp 


I, INTRODUCTION 


Josephson effect has been a fundamental and central 
topic in superconductivity and contributed to determine 
the pairing symmetry in unconventional superconductors 
[lj-[3| . It is well known that standard current-phase rela¬ 
tion (CPR) of Josephson current /(yj) between two super¬ 
conductors is 1(f) ~ sinys, where f is the macroscopic 
phase difference. For d -wave superconductor (DS) junc¬ 
tions, due to the presence of Andreev bound state at the 
interface SHI, exotic quantum interference effects ex¬ 
ist. One is the non-monotonic temperature dependence 
of maximum Josephson current [7-ll| and second is the 
anomalous CPR Due to the presence of ABS, 

sin2y? component of I(ip) is enhanced and free energy 
of the junction can locate neither y? = 0 nor ±7r. Fur¬ 
thermore, a pure sin2y? CPR is possible for d x 2 _ y 2 -wave 
/d xy -wave superconductor j unction [j, El • Thus, d -wave 
junctions have really rich current phase relation and its 
functionalities worth for further research. 

On the other hand, 3D topological insulators (TIs) 
IHl! is a material with a topologically protected sur¬ 
face state due to the strong spin-orbit coupling. The 
generation of superconductivity on the surface state 
of TI via proximity effect has been verified by the 
presence of supercurrent through the Josephson junc¬ 
tions on TI [20|j24j. Josephson current in supercon¬ 
ductor (S)/ferromagnetic insulator(FI)/S junction on TI 
stimulates us since anomalous CPR I(ip) ~ sin(yj — y?o) 
discussed in conventional S/Ferromagnet (F)/S junction 
without TI 25 31] can be realized easily. It is noted that 
fo can be tunable by magnetization [321] . 

We can imagine that more dramatic features will be ex¬ 
pected in d -wave superconductor(DS) /FI/DS Josephson 
junctions on TI. Although there have been several works 
about DS/FI/DS junctions I33l - l35j] . CPR has not been 


clarified for general orientations of the junctions. Recent 
experiments have shown that induced gap function on the 
surface of TI which is formed on high-T c cuprate is almost 
isotropic. This shows that the induced pair potential has 
predominant s-wave symmetry[36]. The possibility of in¬ 
ducing d -wave pairing on TI surface in the actual exper¬ 
iment is still on debate now [37l - l39l ]. Therefore, besides 
the CPR in d -wave superconductor junctions, we must 
study the CPR in junctions with s+d-wave symmetry 
for comparison with actual experiments. 

In this paper, in order to calculate DC Josephson cur¬ 
rent, we develop a formalism of Green’s function of quasi¬ 
particles on the surface of TI jdoj' Both Josephson cur¬ 
rent and local density of states (LDOS) can be calcu¬ 
lated for general orientations of junctions. In general, 
the obtained current phase relation 1(f) has a complex 
if dependence. 1(f) = —I(—f) is easily to be broken 
by magnetization in FI and 1(f) can not be simply ex¬ 
pressed by Iq sin(y? — fo) with nonzero fo- The extreme 
case is d x 2_ y 2 /FI/d xy -wave junctions, where CPR be¬ 
comes sin2y? without magnetization due to the absence 
of the first order Josephson coupling. If we switch on 
magnetization, exotic CPR becomes possible depending 
on the direction of magnetization in FI: i)mixture of 
cosy? and sin(2yj) terms with 1(f) = — I (—f + rd) and 
ii)sin(2yp — 2fo). The complex CPR 1(f) = —I (—f + 7r) 
with 27 t periodicity in i) is not realized in the preexisting 
high-T c cuprate junctions without TI [H-0]. We also cal¬ 
culate Josephson current where s-wave and d-wave pair 
potentials mix. It is found that the anomalous CPR with 
I (f = 0) 0 exists for junctions of asymmetric orienta¬ 

tions even if s-wave component becomes dominant. This 
feature serves as a guide to detect the proximity induced 
d-wave component of pair potential on the surface of TI. 

The outline of this paper is as follows. In section II, we 
present the model and derive the retarded Green’s func- 
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tion and Josephson current. In section III, we show nu¬ 
merical results of Josephson current in DS/FI/DS junc¬ 
tions. In section IV, we show the corresponding results of 
Josephson junctions when s-wave component is induced 
in DS. A concluding remarks will be given in section V. 


II, MODEL AND FORMULAS 



FIG. 1. Schematics of the system: (a) d- wave superconduc¬ 
tor (DS)/ferromagnetic insulator (FI)/DS Josephson junction 
and (b) DS/FI/normal metal (N) junction on the surface of 
a 3D topological insulator. 


As depicted in Fig[TIa), we consider a DS/FI/DS junc¬ 
tion on a 3D TI surface. The effective Hamiltonian for 
the BdG equations is given by 


h(k x ,k y ) + M ia y A( 6 ) 

—io y A* ( 9 ) —h*(—k x , —k y ) — M* ’ 


(1) 


ing region can be expressed as 

Mx) = A ie ik * x + ai A 4 e ik * x + &i A 3 e~ lk * x , ( 2 a) 

Mx) = A 2 e~ ikxX + a 2 A 3 e~ lkxX + b 2 A 4 e ikxX , ( 2 b) 


where k x = fj, cos 9 / Hvf ■ Here the magnitudes of 
the momenta for electrons and holes are approximated 
to be equal since we have made the assumption of 
E, A <C /i. The spinors are given by A 4 = 
(j,e ie ,-e ie 7 i,* 7 i) T , A 2 = (ie i9 ^ 2 , - 72 , 1 , ie i9 ) T , A 3 = 
(ie l9 1 — 1 , 72 , ie l 9 r y 2 ) T and A 4 = (rq, e lf? 7 i,-e* e , i) T 
with 7 t( 2 ) = A 1(2 )/(£’ + \JE 2 - A^ 2) ) and A 1(2) = 

A 0 cos (29 =F 2xi). Other wave functions can be solved 
in a similar way. The retarded green’s function 
G r (x,x',y,y') = J2 ky G l y (*. x ') e ife 43 /-V) can be ob¬ 
tained by combing all the injection processes fioj: 


G r ky (x,x r ) 


' otiipi (x ) (x') + a 2 ipi (x) 'ijjJ (x') + 

a 3 i/j 2 {x) (x') + a 4 ip 2 (x ) (x'), 

{x > x'), 

Pi Ip 3 (X) (x') + p 2 tp 4 (a;) r/>J ( x') + 

P 31 P 3 (®) i> 2 [x!) -(- p 4 i>4 (x) (a/), 

(X < x 1 ), 


( 3 ) 

where ipi— 1^4 are the corresponding conjugated processes 
of ' 0 i = i^ 4 . The coefficients aj = i ^4 and pi = 1^,4 are deter¬ 
mined by satisfying the boundary conditions for all x, x' 
across the regions: 


G r ky {x + 0,x)-G r ky (x-0,x) = h v f (iT z a y ), (4) 


where t x , v , z is the Pauli matrix in the particle-hole space. 
The dc Josephson current for DS/FI/DS junction is de¬ 
termined by electric charge conservation rule 


where h(k x ,k y ) = hvf(k y a x — k x a y ) — /i (0 (—x) + 
0 (x — L)). &x,y,z is the Pauli matrix in the spin 
space and /i is the chemical potential in the su¬ 
perconducting region. The exchange field in FI re¬ 
gion is M = ^i=x,j/,z m *^0 (x) © (L - a?) [s3l - The 
pair potential is given by A 0 (T) cos (29 — 2,\i) 0(— x) + 
A 0 ( T ) cos (29 — 2y 2 ) e~ llp Q(x — A), where ip and 9 are 
the macroscopic superconducting phase and the prop¬ 
agating angle, respectively. The quantity x is taken 
to be the angle between the x-axis and the a-axis 
of d- wave superconductor on top of the TI surface. 
A 0 (T) is assumed to obey the BCS relation A 0 (T) = 
Ao tanh(1.74y / T c /T — 1) with Ao = 1.76fcsT c and T c is 
the critical temperature. 

To construct the retarded Green’s function, we first 
seek the solutions for the four types of quasiparticle injec¬ 
tion processes: left injection for electron (hole): ipi (ip 2 ) 
and right injection for electron (hole): ij ) 3 (0q). Because 
of the translational invariance along y-axis, the wave 
functions 0i( 2 ) = 0 i( 2 ) (x) e tkyV in the left superconduct¬ 


dtP + d x J x + S — 0 , (5) 


where P = e(\E , |T-f- + H/tvp^), j x _ j e x/('I'j’I'| — \P|H/ 7 ) 
and S = 2eIm[A*’F|4'-|- — A*are electric charge 
density, electric current and source term, respectively. 
After straightforward derivation following Ref.[4l[, we 
find that the total Josephson current is given by 
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eksT 
2 h 


^2 sgn(w n )[ 

ky ,UJ n 


Aiai(iuj n ) 

V UJ n + 


A 2 a 2 (iu> n ) 

7 f ; ( 6 ) 

V UJ n + A 2 


< 3 . 1 ( 2 )(iwn) is obtained by analytical continuation E 
to ioj n , where u) n is the Matsubara frequency uj n = 
nksTfin + 1), (n = 0, ±1, ±2....). Eq.© looks simi¬ 
lar to the extended Furusaki-Tsukada’s formulajilj for 
anisotropic d -wave pair potential^, [Hj]. In addition, 
Eq.© is also applicable to the Josephson current of 
s+d wave pairing in which one substitutes A;q 2 ) by 
A 0 + r) A 0 cos (29 2yi) of which 77 > 0 is the ratio be¬ 
tween d -wave pairing and s-wave pairing. 
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III, JOSEPHSON EFFECT IN DS/FI/DS 
JUNCTION 


In this section, we show the results of Josephson cur¬ 
rent I in DS/FI/DS junctions, which has been normal¬ 
ized to eRiyl/Ao where Rn is the interface resistance 
per unit area in the normal state. To analyze the CPR 
further, we decompose the Josephson current into a series 
of different orders of Josephson coupling, 

1 &) = 5Z In sin (nip) + J n cos (nip ), (7) 


where n > 1 is an integer. FigJ^a) shows CPR with¬ 
out magnetization. In this case, the CPR is expressed as 
Y n I n sin (nip) and J n is zero. In the condition with \i = 
0 and X 2 = 7 t/ 4 , the CPR I (ip) becomes Y) n I n sin (nip) 
(n = 2,4...). The feature of this CPR is the same as that 
in the standard d- wave junctions without TI with the pair 
potential considered here. However, as the magnetization 
switches on, the CPR dramatically changes. Figure [^b) 
shows that m y gives a shift of phase difference, which 
is similar to "junctions realized in conventional s-wave 
superconductor/ferromagnet hybrid systems [IB-ED]. As 
the magnetization along m x or ?n 2 -axis appears, the qual¬ 
itative features of CPR of symmetric d x 2_ y 2 /FI/ d x 2_ y 2 
and d xy /Fl/d xy junctions do not change as compared 
to the case without magnetization as shown in dotted 
and dash-dotted lines of Figs. m c) and (d). However, 
in the asymmetric d x 2 _ y 2 /FI/d xy junction, the CPR 
is quite anomalous and the component proportional to 
Y n Jn cos (nip) is generated. We find that I (p) can be 
expressed by Ykl^k sin (2kip) + J 2 fe-i cos( 2 /c—l)y>] where 
k > 1 is an integer, and therefore, I(ip) becomes zero at 
p = ± 7 t /2 as shown in solid lines in Figs. [5Jc) and (d). 
The present CPR is completely different from that of the 
standard d x 2_ y 2 /FI jd xy junction without TI. We can see 
that the term proportional to J\ cos (p) becomes domi¬ 
nant in the limit of large m x or to- in Figs|2](e) and (f). 

To explain the anomalous CPR for nonzero m x or m z 
in d x 2_ y 2/Fl/d xy junction on TI surface, we focus on 
the symmetry of this Hamiltonian. We consider the mir¬ 
ror reflection symmetry with respect to xz- plane, M xz = 
*tJ y To, and the time-reversal symmetry, T = — ia y K,TQ, 
where K. is the complex conjugation operator. In the 
present system, both symmetries are broken. However, 
since the pair potential of d x 2_ y 2 (d xy ) is mirror even 
(odd) with respect to a;z-plane, M xz operation produces 
additional phase, I (p) —» I (p + n). It is also known 
that time reversal operation transforms I (p) to —I (—p)- 
Hence, the composition operator T = M XZ T will give rise 
to I (p) —> —I (—p + ii). Taking into account the fact 
that T makes the (k x ,k y ) state to the (—k x ,k y ) one, we 
can arrive at 


TH(-id x ,k y ,p)f 1 ->• U (id x , k y , -p + n ). ( 8 ) 





cp/jt cp/jt 



FIG. 2. Josephson currents as a function of p in the 
DS/FI/DS junctions. Magnetization in FI is (a) zero, (b) 
along z-axis, (c) along y-axis and (d) along z-axis. Three 
geometries are considered in panels (a)~(d): d x 2_ y 2 /Fl/d xy , 
d x 2 _ y 2 /Fl/d x 2_ y 2 and d xy /Fl/d xy . (e) Josephson currents in 
d x 2_ y 2 /Fl/d xy junctions with m x /p = 0.2, 0.8, and 1.5 and 
(f) those with m z /p = 0.2, 0.8, and 1.5. Other parameters 
are set as T = 0.05T C , p = 1, hvf = 1, A = 0.01 and L = 1. 


It means that —I (—p + tt) = I (p) will be satisfied at any 
p if we consider the junctions between a mirror even and 
mirror odd pair potential. In the d x 2_ y 2/Fl/d xy junction 
with m x or m 2 , we can find that relation (J5J) fulfills at 
any p, which indicates I (p = ± 7 t/ 2 ) = 0 . Above analy¬ 
sis based on mirror reflection symmetry has been applied 
in the Josephson junctions between a singlet and triplet 
superconductor [ 43 . Now, let us look at the Josephson 
current at p = 0 . In the standard DS/FI/DS junc¬ 
tions without TI substrate, due to the spin SU(2) sym¬ 
metry, the rotation or mirror reflection of the ferromag¬ 
netism does not change the CPR and one can always find 
I(p = 0 ) = 00 . However, this SU(2) symmetry is bro¬ 
ken on TI surface due to its nature of spin-momentum 
locking and thus I(p = 0 ) becomes nonzero which gen¬ 
erates exotic 27r-periodic CPR —I(—p + n) = I(p). 

Next, we plot the temperature dependence of the 
maximum Josephson current I c of DS/FI/DS junc¬ 
tions in the left panels of FigUl For simplicity, 
only the z component of magnetization m z is consid¬ 
ered. We concentrate on the low temperature region 
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FIG. 3 . (a) The maximum Josephson currents in the 

d x 2_ y 2/Fl/d x 2_ y 2 junctions, (b) The LDOS on the surface 
of DS in the d x 2_ y 2 /Fl/N junctions, po is the electron den¬ 
sity of states of the bulk N at Fermi energy, (c) and (d) are 
the maximum Josephson currents and LDOS in d xy /Fl/d xy 
junctions, respectively. The magnetization is along a-axis in 
all panels. Other parameters are set as the same as Fig( 2 ] 


T/T c < 0.4 in which the behavior of I c is highly 
influenced by the zero-energy states (ZESs). There¬ 
fore, we display the surface density of states at the 
edge of DS. It is obtained by calculating the LDOS 
p e {x,E) = + G r ky 22 (x, x, E)\ 

at the DS/FI interface in the DS/FI/N junction as il¬ 
lustrated in FigQJb). From Fig[3[a), we can see that 
temperature dependence of I c in d x 2 _ y 2 /FI/d a , 2_, [ 2 junc¬ 
tions changes from Kulik-Omelyanchuk (K-0)[45[ type 
to Ambegaokar-Baratoff (A-B)[46j type with decreasing 
m z . However, in d xy /Fl/d xy junction, this tendency is 
reversed when we decrease to,, as shown in Fig[3tc). As 
seen from Figs|3|b) and (d), we can see that as LDOS at 
zero energy is enhanced, temperature dependence of I c 
is reduced to be the K-0 type. In other case, it is in the 
A-B type. Since the spin degeneracy is lifted on the sur¬ 
face states of TI by spin momentum locking, we obtain 
the hig hly asymmetric Yu-Shiba-Rusinov type of LDOS 
47 -491 in the d r 2_ v 2 /FI interface. This finding is similar 
to that in the s-wave superconductor/FI/N j unction [50j. 


IV, JOSEPHSON EFFECT WITH s+d-WAVE 
PAIRING 


Recent experiments have shown that the induced en¬ 
ergy gap by high T c cuprate on the surface of TI is al¬ 


most isotropic 


It is interesting to clarify the role 


of the induced s-wave pair potential on DS/FI/DS junc¬ 
tions. In this section, we calculate Josephson current in 


s+d-wave/FI/s+d-wavejunctions on TI. The pair poten¬ 
tial is A = A 0 + 77 A 0 cos(20 — 2\) with \ = 0(7t/ 4) on 
the left(right) side. The ratio 77 is chosen to be 0.5 so 
that the system is s-wave dominant and fully gapped. 
The obtained I(ip) has a typical sinusoidal shape of s- 
wave Josephson current where the first order coupling 
I\ sin ip plays the predominant role. Because of the s- 
wave component of pair potential, the mirror reflection 
symmetry T at ip = ±7t/ 2 is broken and thus nonzero 
current I(ip = ±7t/ 2) can be expected. Also, in the pres¬ 
ence of m x or m z , we find a nonzero Josephson current 
at ip = 0 , 7r in such junctions. The obtained anoma¬ 
lous CPR in the s+d-wave Josephson junctions can be 
used to probe the d -wave component of the induced pair 
potential on TI surface. For example, one can observe 
the supercurrent flow without macroscopic phase differ¬ 
ence in s+d x 2 _ y 2 /FI/ s+d xy junctions. As seen in FigjTj 
the existence of d -wave component generates a nonzero 
current I(p = 0) when one turns on either m x or m z . 



FIG. 4 . I(pp = 0 ) as a function of 77 in s+d x 2_ y 2/FI/s+d xy 
junctions on TI. Other parameters are set as the same as in 
FigEJ 


V, CONCLUSION 

In summary, we have theoretically studied the Joseph¬ 
son effect in d -wave superconductor-ferromagnet insula¬ 
tor (FI) hybrids on the surface of TI. Depending on the 
orientation of the magnetization in FI, the exotic current- 
phase relation which violates I (ip) 7 ^ —I (—ip) have been 
obtained in two different ways: (i) through a simple phase 
shift and (ii) mixture of cos ip term into the original CPR. 
The latter case can generate the exotic current-phase 
relation I(ip) = —I(—p + 7 r) with 27r periodicity. We 
show that the Josephson current is enhanced due to the 
zero energy states on the edge of d -wave superconduc¬ 
tor. For comparison with actual experiments, we calcu¬ 
late the Josephson current when both s- and d -wave pair 
potentials exist. The anomalous current-phase relation 
is also found which provides a way to probe the finger¬ 
print of d -wave pair potential in high-T/ cuprate/TI het¬ 
erostructures. Our preliminary theoretical investigation 
has practical significance for controlling the Josephson 
current and designing new functional devices. 
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